Charge transport in underdoped bilayer cuprates 
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Within the t-J model, we study the charge transport in underdoped bilayer cuprates by considering 
the bilayer interaction. Although the bilayer interaction leads to the band splitting in the electronic 
structure, the qualitative behavior ol the charge transport is the same as in the case of single layer 
cuprates. The conductivity spectrum shows a low-energy peak and the unusual midinfrared band. 
This midinfrared band is suppressed severely with increasing temperatures, while the resistivity in 
the heavily underdoped regime is characterized by a crossover from the high temperature metallic- 
like to the low temperature insulating-like behaviors, which are consistent with the experiments. 
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It has become clear in the past ten years that cuprate 
superconductors are among the most complex systems 
studied in condensed matter physics 1,2 . The compli- 
cations arise mainly from (1) strong anisotropy in the 
properties parallel and perpendicular to the Cu02 planes 
which arc the key structural element in the whole cuprate 
superconductors, and (2) extreme sensitivity of the prop- 
erties to the compositions (stoichiometry) which control 
the carrier density in the CUO2 plane, and therefore the 
regimes have been classified into the underdoped, opti- 
mally doped, and overdoped, respectively 1 ' 2 . In the un- 
derdoped and optimally doped regimes, the experimental 
results 3 show that the ratio of the c-axis and in-plane re- 
sistivities R — p c (T)/p a b(T) ranges from R ~ 100 to 
R > 10 5 , which reflects that the charged carriers are 
tightly confined to the Cu02 planes. This large magni- 
tude of the resistivity anisotropy also leads to the gen- 
eral notion that the physics of doped cuprates is almost 
entirely two-dimensional, and can be well described by a 
single Cu02 plane 4 . However, this picture seems to be in- 
compatible with the fact that the superconducting tran- 
sition temperature T c is closely related to the number of 
CuC>2 planes per unit cell, with single layer compounds 
of a family generically having lower T c than bilayer or 
trilayer compounds 2 . Additionally, there are some sub- 
tle differences of the magnetic behaviors between doped 
single layer and bilayer cuprates. By virtue of system- 
atic studies using NMR and pSR techniques, particu- 
larly the inelastic neutron scattering, only incommensu- 
rate neutron scattering peaks for the single layer lan- 
thanum cuprate are observed in the underdoped regime 5 , 
however, both low-energy incommensurate neutron scat- 
tering peaks and high-energy commensurate [7r,7r] reso- 



nance for the bilayer yttrium cuprate in the normal state 
are detected 6 . These experimental results highlight the 
importance of some sort of coupling between the Cu02 
planes within a unit cell. It is believed that all these 
experiments produce interesting data that introduce the 
important constraints on the microscopic models and the- 
ories. 

The charge transport of doped single layer cuprates has 
been addressed from several theoretical viewpoints 7,8 . 
Based on the charge-spin separation, an attractive pro- 
posal is spinons and holons as basic low-energy excita- 
tions, serving as the starting point for the gauge-theory 
approach 7 . It has been shown 7 within the t-J model that 
above the Bose-Einstein temperature, the boson inverse 
lifetime due to scattering by the gauge field is of order 
T, which suppresses the condensation temperature and 
leads to a linear T resistivity. On the other hand, the 
spin-fermion model near the antiferromagnetic instability 
has been developed to study the normal-state properties 
of doped cuprates 8 . This spin-fermion model describes 
low-energy fermions interacting with their own collective 
spin fluctuations. Within this approach 8 , the anomalous 
transport of doped single layer cuprates has been stud- 
ied extensively 8 , and the results are consistent with the 
experiments. 

As regards an intracell hopping, the band splitting 
in doped bilayer cuprates was shown by the band 
calculation 9 , and clearly observed 10,11 recently by the 
angle-resolved-photoemission spectroscopy in the doped 
bilayer cuprate BiaSraCaCuaOg+a above T c . This bilayer 
band splitting is due to a nonvanishing intracell coupling. 
Moreover, the magnitude of the bilayer splitting is con- 
stant over a large range of dopings 11 . Considering these 



1 



highly unusual normal state properties in the underdoped 
regime 1 ' 2,5 ' 6 , a natural question is what is the effect of 
the intracell coupling on the normal state properties of 
doped bilayer cuprates. This is a challenge issue since the 
mechanism for the superconductivity in doped cuprates 
has been widely recognized to be closely related with the 
anisotropic normal-state properties 12 . Based on the t-J 
model, the charge transport and spin response of doped 
single layer cuprates in the underdoped regime have been 
discussed 13-15 within the fermion-spin theory 16 , and the 
obtained results are consistent with experiments 17 . In 
this paper, we apply this successful approach to study 
the charge transport of the underdoped bilayer cuprates. 
Our results show that although the bilayer interaction 
leads to the band splitting in the electronic structure, 
the qualitative behavior of the conductivity and resistiv- 
ity is the same as in the single layer case. The conduc- 
tivity shows the non-Drude behavior at low energies and 
anomalous midinfrared band separated by the charge- 
transfer gap, while the temperature dependent resistivity 
in the heavily underdoped regime is characterized by a 
crossover from the high temperature metallic-like to the 
low temperature insulating-like behaviors. 

We start from the bilayer t-J model, which can be 
written as, 

H = - f E - *l EccLcw + h.c.) 

aii)<7 icr 
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With J cff = J[(l - 5) 2 - 4>% J±eff = J±[(l - 5f - (ft], 

the holon particle-hole order parameters <j> = (/i^/i a j+^) 
and cf>± = (h\ i Ii2i), S is the hole doping concentra- 
tion, and S^i (Sai) is the pseudospin raising (lowering) 
operator. Since the single occupancy local constraint 
has been treated properly within the fermion-spin the- 
ory, this leads to disappearing of the extra gauge de- 
gree of freedom related with this local constraint under 
the charge-spin separation 16 . In this case, the charge 
fluctuation couples only to holons 13 ' 14 . However, the 
strong correlation between holons and spinons is still in- 
cluded self-consistcntly through the spinon's order pa- 
rameters entering the holon's propagator, therefore both 
holons and spinons are responsible for the charge trans- 
port. In this case, the conductivity can be expressed 
as cr(u>) = -Imll^Ww, with nW(a;) is the holon 
current-current correlation function, and is defined as 
n(>0( T _ T ') = -(T T j( h \T)j( h \T')), where r and r' are 
the imaginary times, and T T is the r order operator. 
Within the Hamiltonian (2), the current density of holons 
is obtained by the time derivation of the polarization op- 
erator using Heisenberg's equation of motion as, 
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where fj = ±x, ±y within the plane, a = 1 and 2 is 
plane indices, C\ ia (C a i a ) is the electron creation (an- 
nihilation) operator, S ai = are spin opera- 
tors with a = (a x ,(jy,a z ) as Pauli matrices, and fi is 
the chemical potential. The bilayer t-J model (1) is 
defined in the subspace with no doubly occupied sites, 
*- e -' S<7 C aia C a i a < 1. The strong electron correla- 
tion in the t-J model manifests itself by this single oc- 
cupancy on-site local constraint 4 . To deal with the lo- 
cal constraint in analytical calculations, the fermion-spin 
theory 16 , C a% \ = h^S^ and C m \ = h^S^, has been 
proposed, where the spinless fermion operator h ai keeps 
track of the charge (holon), while the pseudospin opera- 
tor S a i keeps track of the spin (spinon), then it naturally 
incorporates the physics of the charge-spin separation. 
In this case, the low-energy behavior of the bilayer t- 
J model (1) in the fermion-spin representation can be 
rewritten as, 

H = 1 E h li+fi hal ( S ai S ai+fi + S ai S li+f)) 
aii) 

+ t ± 5>U 2 > + hihu^S+S^ + S U S+) 



where Ru (i?2i) is lattice site of the Cu02 plane 1 (plane 
2), X = (Sai S ai+fi) and X± = (SuSx) are the spinon 
correlation functions, and e is the electronic charge, 
which is set as the unit hereafter. The holon current- 
current correlation function can be calculated in terms 
of the holon Green's function g(k,uj) as in the single 
layer case 13,14 . However, in the bilayer system, be- 
cause there are two coupled Cu02 planes, then the en- 
ergy spectrum has two branches. In this case, the one- 
particle holon Green's function can be expressed as a 
matrix g(i - j,r - t') = g L (i — j,r - t') + a x g T {i - 
j, t — t'), with the longitudinal and transverse parts are 
defined as g L (i - j,r - r') = -(T T /i ai (T)/i^-(r')) and 
gr(i - j,r- t') = -{Trhai^hl^ir')) {a ^ a'), respec- 
tively. Following discussions of the single layer case 13 ' 14 , 
we obtain the conductivity of doped bilayer cuprates as 
cr(u>) = a^ L \u>) + <t( t )(w), with the longitudinal and 
transverse parts are given by, 



^AW{k,u' + u>)A<£\k,u')x 



2 



n F (uj' + u) — np{io') 



(4a) 



k 



\{2Z X t lsk ) 2 - (2 X± t ± ) 2 } x 



■/-oo 27r 

njr(w' + w) — % (w') 



(4b) 



respectively, where Z is the coordination number within 
the plane, 7 s fc = (sinfcj, + sinfc a )/2, and uf(lj) is 
the fermion distribution function. The longitudi- 
nal and transverse holon spectral functions are ob- 
tained as Aj/{k,w) = —2lmgL(k,uj) and A^\k,w) = 
— 2lmgT(k, u>), respectively. The full holon Green's 
function g^^k,^) = g ( -^- 1 (k,iv) - Y,( h \k,u) with 
the longitudinal and transverse mean-field (MF) holon 

Green's functions, g^\k,u>) — 1/2 l/( w — $f*) an d 

$\k,u) = l/2E,(-l)" +1 /(w - tP), where v = 1, 2, 
and the longitudinal and transverse second-order holon 
self-energy from the spinon pair bubble are obtained by 
the loop expansion to the second-order as, 
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respectively, with E„„'„"(fc,p, q, uj) is given by, 
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with 7 fe = (cosfc x + cosfc,,)/2, B^- 1 = B k - Jxeff[X-L + 
2xl(-l)1[e± + (-in Sfe - A[(2e X *+ X ) 7fe -(e X +2 X *)], 
A = 2ZJ eff , e = 1 + 2t(j)/J cS , ej_ = 1 + 4tj_<j>j_/ Jj_ e ff, and 

+ n B ( W J2)[l + n B (o;M)] ) 
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~ M&'2)[1 + + » B (f)], (7) 

where tib(wJ^) is the boson distribution function, the 

MF holon excitation ^ = 2Z X tj k + n + 2 Xl _t ± (-l) v+1 , 

the MF spinon excitation (w^) 2 = w 2 . + A|(-l)" +1 , 
with = A^l + A 2lk + A 3 , A 2 = X^ k + X 2 , and 

A 1 = ae\ 2 (x/2 + e X z ), 

A 2 = e\ 2 {(l-Z)a(e X /2 + X z )/Z 

- a(C z + C/2)- (l-a)/(2Z)] 

- aAJ ±off [e(Ci + xl) + e±(<?j_ + eX±)/2], 
A 3 = \ 2 [a(C z + e 2 C/2) + (1 - a)(l + e 2 )/(AZ) 

- ae( X /2 + e X z )/Z] 

+ a\J ±cS [ee ± C ± + 2C Z ± ] + J 2 ±cB (el + l)/4, 
X 1 = aAJ ±off [(e ± x + e X± )/2 + ee ± ( X z ± + X z )], 
X 2 = -a\J ±cS [ee ±X /2 + e ± ( X z + Cj_) + eC ± /2] 

- ^XeffA (8) 

where the spinon correlation functions X z — (S^S^-), 
Xl = (S z u S z 2i ), C = {l/Z 2 )Y,ri>K i+n S-^,), and 
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a/Z)^^-^), and C'l = (l/Z^iS'At). In 
order to satisfy the sum rule for the correlation function 
(S'j" i S , ~) = 1/2 in the absence of the antiferromagnetic 
long range order, a decoupling parameter a has been in- 
troduced in the MF calculation, which can be regarded as 
the vertex correction 18 . All these order parameters, de- 
coupling parameter a, and the chemical potential /i have 
been determined self-consistently, as done in the single 
layer case 18 . 

The frequency- and temperature-dependent conductiv- 
ity is a powerful probe for systems of interacting elec- 
trons, and provides very detailed informations of the exci- 
tations, which interacts with carriers in the normal-state 
and might play an important role in the superconduc- 
tivity. In Fig. 1, we present the results of the conduc- 
tivity a(oj) at doping S = 0.05 (solid line), S = 0.06 
(dashed line), and S — 0.07 (dotted line) for parameters 
t/J = 2.5, tx/t = 0.25, Jl/J = 0.25 with tempera- 
ture T = in comparison with the experimental data 19 
taken on the underdoped YBa2Cu307-a; (YBCO) (inset). 
The conductivity of bilaycr cuprates in the underdoped 
regime shows a sharp low-energy peak at uu < 0.5£ and the 
unusual midinfrarcd band appearing inside the charge- 
transfer gap of the undoped system. After an analy- 
sis, we found that this low-energy peak decays fastly as 
cr(w) ~ 1/co (non-Drude fall-off) with increasing energies. 
Moreover, the weight of the midinfrared peak is doping 
dependent, and the peak position is shifted to low energy 
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with increasing dopings. For a better understanding of 
the optical properties of doped bilayer cuprates, we have 
studied conductivity at different temperatures, and the 
results at doping S = 0.06 for t/J = 2.5, t±/t = 0.25, 
and J ± /J = 0.25 in T = (solid line), T = 0.3J (dashed 
line), and T = 0.5J (dotted line) are plotted in Fig. 2 
in comparison with the experimental data 19 taken on 
the underdoped YBCO (inset). It is shown that a(uj) 
is temperature dependent, and the charge-transfer gap 
is severely suppressed with increasing temperatures, and 
vanishes at higher temperature (T > 0.4J). Our results 
are in qualitative agreement with the experiments 19 . In 
comparison with the results from Refs. 13,14 , it is shown 
that the present conductivity also is qualitatively con- 
sistent with these in the single layer case. In the above 
calculations, we also find that the conductivity a(u>) is es- 
sentially determined by its longitudinal part <j( l \uj), this 
is why in the present doped bilayer cuprates the conduc- 
tivity spectrum appears to reflect the single layer nature 
of the electronic state 1 ' 2,13,14 . This is also why the in- 
plane charge dynamics is rather universal within whole 
doped cuprates 1 ' 2 . 

Now we turn to discuss the resistivity, which is closely 
related to the conductivity, and can be obtained as 
p(T) = l/lim^^o c(w). This resistivity has been cal- 
culated, and the results at doping 5 = 0.05 (solid line), 
S = 0.06 (dashed line), and S = 0.07 (dotted line) for 
t/J = 2.5, t±/t = 0.25, and J ± /J = 0.25 arc plotted 
in Fig. 3 in comparison with the experimental results 20 
taken on the underdoped YBCO (inset). These results 
show that in the heavily underdoped regime, although 
the temperature-dependent resistivity is characterized by 
a crossover from the high temperature metallic-like to 
the low temperature insulating-like behaviors, the nearly 
temperature liner dependence in the resistivity domi- 
nates over a wide temperature range, in agreement with 
the experimental results 20 . In comparison with the re- 
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FIG. 1. The conductivity at <5 = 0.05 (solid line), 8 = 0.06 
(dashed line), and 5 = 0.07 (dotted line) for t/J = 2.5, 
t±/t = 0.25, and J±/J = 0.25 in the zero temperature. Inset: 
the experimental result on the underdoped YBa2Cu307_ 2! 
taken from Ref. 19 . 



suits from Refs. 13 ' , it is shown that the present resistiv- 
ity also is qualitatively consistent with these in the single 
layer case. We emphasize that since the order parame- 
ters, decoupling parameter a, and the chemical poten- 
tial [i have been determined self-consistently, then these 
theoretical results were obtained without any adjustable 
parameters. Furthermore, it is found in the above dis- 
cussions that the present results are insensitive to the 
reasonable values of t/J, t±/t, and J±/J as in the single 
layer case 13 ' 14 . 

An explanation for the metal-to-insulating crossover in 
the resistivity in the heavily underdoped regime can be 
found from the competition between the kinetic energy 
and magnetic energy in the system. Since cuprate su- 
perconducting materials are doped Mott insulators, ob- 
tained by chemically adding charge carriers to a strongly 
correlated antiferromagnctic insulating state, therefore 
doped cuprates are characterized by the competition be- 
tween the kinetic energy (t) and magnetic energy (J). 




FIG. 2. The conductivity at 8 = 0.06 for t/J = 2.5, 
t±/t = 0.25, and J ± /J = 0.25 in T = (solid line), T = 0.3J 
(dashed line), and T = 0.5J (dotted line). Inset: the exper- 
imental result on the underdoped YBa2Cu307_a; taken from 
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FIG. 3. The resistivity at 5 = 0.05 (solid line), 8 = 0.06 
(dashed line), and 8 = 0.07 (dotted line) for t/J = 2.5, 
t±/t — 0.25, and J±/J = 0.25. Inset: the experimental result 
on the underdoped YBa2Cu307_i, taken from Ref. 20 . 
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The magnetic energy J favors the magnetic order for 
spins, while the kinetic energy t favors derealization of 
holes and tends to destroy the magnetic order. In the 
present fermion-spin theory, although both holons and 
spinons contribute to the charge transport, the scatter- 
ing of holons dominates the charge transport 13 , where 
the charged holon scattering rate is obtained from the 
full holon Green's function (then the holon self-energy 
(5) and holon spectral function) by considering the holon- 
spinon interaction, therefore in the heavily underdoped 
regime the observed crossover from the high temperature 
metallic-like to the low temperature insulating-like be- 
haviors in the resistivity is closely related with this com- 
petition. In lower temperatures, the holon kinetic energy 
is much smaller than the magnetic energy, in this case 
the magnetic fluctuation is strong enough to severely re- 
duce the charged holon scattering and thus is responsible 
for the insulating-like behavior in the resistivity. With 
increasing temperatures, the holon kinetic energy is in- 
creased, while the spinon magnetic energy is decreased. 
In the region where the holon kinetic energy is much 
larger than the spinon magnetic energy at higher tem- 
peratures, the charged holon scattering would give rise 
to the temperature linear resistivity. 

In summary, we have studied the charge transport in 
the underdoped bilayer cuprates by considering the bi- 
layer interaction. It is shown that although the bilayer 
interaction leads to the band splitting in the electronic 
structure, the qualitative behavior of the charge trans- 
port is the same as in the single layer case. The conduc- 
tivity spectrum shows a low-energy peak and the anoma- 
lous midinfrared band. This midinfrared band is sup- 
pressed severely with increasing temperatures, while the 
resistivity exhibits a crossover from the high temperature 
metallic-like to the low temperature insulating-like be- 
haviors. Our results also show that the mechanism that 
cause this unusual charge transport in the underdoped 
cuprates is closely related to the background antiferro- 
magnetic correlations. 
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